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Abstract 

We apply lattice Boltzmann methods to study the relaxation of the velocity 
profile in binary fluids under shear during spinodal decomposition. In simple 
fluids, when a shear flow is applied on the boundaries of the system, the time 
required to obtain a triangular profile is inversely proportional to the viscosity 
and proportional to the square of the size of the system. We find that the 
same behaviour also occurs for binary mixtures, for any component ratio in 
the mixture and independently from the time when shear flow is switched on 
during phase separation. 
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I. INTRODUCTION 



Recently, the dynamical behaviour of complex fluid and binary mixtures under the action 
of applied flows has been investigated in many theoretical and experimental works [1]. In 
phase separating binary mixtures the morphology and the growth properties of the domains 
of the two components are greatly affected by the presence of a flow [2] . This flow is generally 
imposed from the boundaries onto the system and takes some time to reach its steady value. 
If this time is comparable with other typical time-scales in the system, transient effects have 
to be carefully considered. For example, in the case of an oscillating shear flow, the ratio 
between the period of the applied flow and the relaxation time of a triangular steady shear 
profile is relevant for the effectiveness of the flow on the bulk of the system. In spinodal 
decomposition, if that ratio is small, the growth remains mainly isotropic and elongated 
domains can be observed only in the boundary layers [3]. 

The estimation of relaxation times of applied flows is generally based on analogies with 
simple fluids. However, a careful checking of the laws valid for simple fluids in cases when 
mesoscopic structures are present in the system - e.g. interfaces between domains in spinodal 
decomposition, cannot be found in literature. In this paper we evaluate numerically the re- 
laxational properties of a triangular shear flow applied on phase separating binary mixtures. 
At a certain time during phase separation the top and the bottom boundaries of our system 
start to move with velocities w and — w respectively and after a while a triangular profile 
will appear. We study the evolution of the velocity field in terms of the viscosity and of the 
size of the system and find that the results for a simple fluid also hold in the case of binary 
mixtures with interfaces inside. 

We simulate the binary mixture system by applying lattice Boltzmann methods which 
allow to mimic the behaviour of a system described by Navier-Stokes and convection-diffusion 
equations. The lattice Boltzmann algorithm is based on a collision step and a propagation 
step occurring on time At. Only a finite sets of local velocities are allowed. In particular we 
use a scheme based on a free-energy approach which has the advantage that the equilibrium 
thermodynamics of the system is "a priori' known. The role of the simulation time step At 
will be also considered. With fast relaxing velocity profiles, the usual choice At = 1 is not 
appropriate for obtaining smoothly relaxing profiles and smaller At have to be considered. 

In the next section we briefly describe the methods used in the simulations; section III 
contains our results and section IV some conclusions. 

II. THE MODEL 

We consider a two-dimensional binary fluid with components A and B of number den- 
sity pa and respectively. Such a system can be modeled by the following free energy 
functional, 

F = Jd r [ 1 -p\np+^+ b -^ + ^n (1) 

where p = Pa+Pb is the local total density and (p = pa~Pb is the local density difference and 
the order parameter of the system; the term in p gives rise to a positive background pressure 
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and does not affect the phase behavior. The terms in ip correspond to the usual Ginzburg- 
Landau free energy typically used in studies of phase separation [4]. The polynomial terms 
are related to the bulk properties of the fluid. The gradient term is related to the interfacial 
properties. The parameter b is always positive, while the sign of a distinguishes between 
a disordered (a > 0) and a segregated mixture (a < 0) where the two pure phases with 

p = ±^J—a/b coexist. In this paper we will consider quenches into the coexistence region 
with a < and b = —a, so the equilibrium values for the order parameter are <p — ±1. 
The initial state in simulations will be random configurations corresponding to the high 
temperature disordered phase. 

The thermodynamic properties of the fluid follow directly from the free energy (1). The 
chemical potential difference between the two fluids is given by 

5F 

Afi = — = ap + bp 3 - reVV (2) 
op 

The pressure is a tensor P a p since interfaces in the fluids can exert non-isotropic forces [6]. 
A suitable choice is 

Pa(3 = PoSap + Kdapdpip, (3) 

where the diagonal part po can be calculated using thermodynamics relations from (1): 

5F 5F 

= p - + a -P 2 + ^-^p)-^iypf (4) 



Our simulations are based on the lattice Boltzmann scheme developed by Orlandini et 
al [8]. and Swift et al. [9] . In this scheme the equilibrium properties of the system can 
be controlled by introducing a free energy which enters properly into the lattice Boltzmann 
model. The scheme used in this paper is based on the D2Q9 lattice: A square lattice is used 
in which each site is connected with nearest and next-nearest neighbors. The horizontal and 
vertical links have length Ax, the diagonal links \^2Ax, where Ax is the space step. Two sets 
of distribution functions fi(r,t) and gi(r,t) are defined on each lattice site r at each time t. 
Each of them is associated with a velocity vector e^. Defined At as the simulation time step, 
the quantities ejAt are constrained to be lattice vectors so that |ej| = Ax/ At = c for % — 1, 
2, 3, 4 and |ej| = \[2c for i — 5, 6, 7, 8. Two functions fo(r,t) and g (r,t), corresponding 
to the distribution components that do not propagate (e = 0), are also taken into account. 
The distribution functions evolve during the time step At according to a single relaxation 
time Boltzmann equation [10,11]: 

fi(r + e t At, t + At)- f t (r, t) = --[/,(r, t) - /f(r, t)], (5) 



&(r + e,At, t + At)- 9i (r, t) = [^(r, t) - <?f (r, *)], (6) 
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where r and t 9 are independent relaxation parameters, f^ 9 (r,t) and g^ q (r,t) are local equi- 
librium distribution functions. Following the standard lattice Boltzmann prescription, the 
local equilibrium distribution functions can be expressed as an expansion at the second order 
in the velocity v [12,13]: 

f? = A + C v 2 

f- q = A r + Biv a e ioi + Cjv 2 + D I v a vpe ia e i p + G Iia pe ia e i/3 
i = l,2,3,4, 

(7) 

f- g = An + B n v a e ia + C n v 2 + D n v a vpe ia eif3 + G n 
i = 5, 6, 7, 8, 

and similarly for the g e9 , % = 0, 8. The expansion coefficients A Q , Ai, An, Bi, ■ ■ ■ are 
determined by using the following relations 

E*/fM)=p, E^f(r, i) = (8) 
Ei fi 9 ( r , t)et = pv, Ei ^ e9 e ia = ipv a , (9) 

Ei fi q e ia e^ = c 2 P afS + pv a V{3, Ei 9?e-ia.e-i{i = c 2 TA/j,5 af s + (pv a Vf3, (10) 

where P aj « is the pressure tensor, A/x is the chemical potential difference between the two 
fluids and T is a coefficient related to the mobility of the fluid. We stress that we are con- 
sidering a mixture with two fluids having the same mechanical properties and, in particular, 
the same viscosity. The second constraint in Eq. (9) expresses the fact that the two fluids 
have the same velocity. 

A suitable choice of the coefficients in the expansions (7) is shown in Ref. [3]. Such 
a lattice Boltzmann scheme simulates at second order in At the continuity, the quasi- 
incompressible Navier-Stokes and the convection-diffusion equations with the kinematic vis- 
cosity v and the macroscopic mobility given by [12-14] 

z/ = Atf(r-i), Q = TAtc 2 (r^-l). (11) 

The shear flow can be imposed by introducing boundary walls on the top and bottom 
rows of lattice sites. The velocities of the upper and lower walls are along the horizontal 
direction and their values are w and —w, respectively, where w = ^{L — l)/2 and 7 is 
the shear rate imposed on the system. The bounce-back rule [15,16] is adopted for the 
distribution functions normal to the boundary walls. In order to preserve correctly mass 
conservation, a further constraint, related to the distribution functions at the previous time 
step (t — At), is used. Details of the scheme are given in Ref. [17]. 

Finally, we observe that if we set a = b = k = in the free energy functional (1), the 
present lattice Boltzmann methods can be used to simulate simple fluids. 
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III. RESULTS 

We numerically check the relaxation behaviors of the horizontal velocity profile for sym- 
metric and asymmetric binary fluids. We consider two cases: (i) switch on the shear from 
the beginning of the phase separation and (ii) switch on the shear during the phase sepa- 
rating process and after the interfaces between the two fluids have formed. We denote the 
time at which the shear is switched on as t on . We focus on the effects of viscosity, so we 
vary r in the simulations. The other parameters, if not differently stated, are fixed at the 
values: -a = b = 1.252 x 10" 4 , k = 8 x 10" 5 , T = 40, Ax = 1, At = 0.2, 7 = 0.005. Similar 
results are obtained for other sets of parameters. 

Figure 1 shows the relaxation of the horizontal velocity profile for a symmetric binary 
fluid with r = 2, L = 64, t on = 5000. The four lines correspond to t = 5020, t = 5040, 
t = 5060, t = 5120. From this figure we can see how the shear flow comes into the bulk of the 
system. When t > 5120 the velocity profile shows an almost linear behavior with y. Before 
switching on the shear the domains grow isotropically. After t on anisotropic behavior come 
into the bulk of the system with time. Figure 2 shows the configuration at the time t = 5120. 
Domains separated by well formed interfaces can be observed to incline to the flow direction 
with time. For the case t on — we observed the same behavior for the velocity profile, while 
the interfaces have not reached their equilibrium shape. That means the existence of the 
interfaces does not influence the shear effects coming into the system. For the asymmetric 
case, we studied the binary fluid with the ratio 80 : 20. After switching on the shear the 
velocity profile shows the same behavior. 

To understand better the time evolution of the velocity profile of binary mixtures, we 
take the simple fluid as a reference. We consider the Newtonian viscous flow between two 
infinite plates with a distance of L. We use y to denote the coordinate in the vertical 
direction and y = in the middle of the system. The velocities of the upper and the lower 
plates are w and — w, respectively . The shear rate imposed on the system is 7. The motion 
equation is 

du d 2 u 
P dt =V dy^ 

The velocity profile can be obtained by standard methods [18] and is given by 

u = iy- EC" 1 ) + — ex P( Fa— *) M-rV)- (13) 

n H7T Li Li 

When t is large enough, the modes with n > 2 can be neglected, which is confirmed by 
simulations. So we can define a relaxation time for the velocity profile in the following way, 

T - L " - 
R 4ir 2 v 27r 2 Atc 2 (2r- 1) 

It is interesting to see if or not such a definition also works for binary fluids. To numerically 
check the relation between Tr and v (or r) or L, we calculate Tr in the following way, 



(12) 



(14) 



where 



^y~ u ) ( i 6) 

In order to calculate Tr, we use a time t at which the velocity profile has been almost linear, 
so that terms with n > 2 in Eq.(13) can be neglected. 

Figure 3 shows the simulation results for Tr as a function of u, where L = 128, At = 0.2, 
the vertical axis is for 1/Tr and the horizontal axis is for v so that we expect a linear 
behavior. The dotted line with points correspond to the simulation results and the solid line 
corresponds to the expected value from the definition (14). The simulation results confirm 
the existence of the exponential behavior in the relaxation process of the velocity profile and 
the validity of the definition of Tr. 

Figure 4 shows Tr as a function of L, where the vertical axis is for Tr and the horizontal 
axis is for L? . The dotted lines with symbols are simulation results for the cases of v — 4.17, 
v = 5.83 and v = 7.5. The solid lines are the expected values from the definition (14). The 
four points in each case correspond to L — 64, 128, 256, and 512. The simulation step in 
this figure is At = 0.2. We find the expected linear behavior between Tr and L? . 

There is a bending tendency in the simulation results for 1/Tr in Fig. 3, so that the 
simulation results deviate from the expected values. If we continue the simulation up to 
a much higher viscosity region, we will find a more pronounced bending behavior. We 
emphasize that this is an artificial phenomenon, which can be mainly attributed to the 
finite size of the simulation step At. If the viscosity is larger the velocity profile relaxes more 
quickly and we should use a smaller time step At for the simulation. When the viscosity 
is very high and At is not small enough we can not observe smooth velocity profiles in the 
simulation. To confirm this numerical analysis we show two simulation results in Fig.5, 
where the definition value of 1/Tr is also shown to guide the eyes. Compared with the case 
of At = 0.2, when we use At = 0.1 the linear behavior region is almost doubled. 

Finally, it is interesting to compare the relaxation behaviors of binary and simple fluids. 
Fig. 6 shows the simulation results of the velocity profiles for the case with L = 256, r = 2 
(v = 2.5), t — t on = 600. The solid line corresponds to Eq. (13) with n — \. From this figure 
the following remarks are evident: (i) the relaxation process of the binary mixtures follows 
the same behavior as that of simple fluids; (ii)for binary mixtures the relaxation behavior 
of the velocity profile is independent of the component ratio and the time when shear is 
switched on, which means that the forming of the interfaces does not evidently influence the 
shear effects coming into the system. 



IV. CONCLUSIONS 

In this paper we study the relaxation behavior of the velocity profile of binary mixtures 
under steady shear with lattice Boltzmann methods. Following the simple Newtonian vis- 
cous flow, we define a relaxation time of the velocity profile Tr whose value is inversely 
proportional to the viscosity and proportional to the square of the size of the system. The 
simulation results show that the shear effects come into the system in the same way for 
binary mixtures and simple fluids, which confirm the validity of the definition of Tr in this 
and previous studies [3]. For binary mixtures, the shear behavior is independent of the 
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component ratio and the time at which the shear is switched on. The presence of inter- 
faces between the two fluids has negligible influence on the relaxation process of the velocity 
profile. 
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FIG. 1. Relaxation process of the horizontal velocity profile for the symmetric binary fluids, 
where r = 2, L = 64, t on = 5000. The corresponding times are shown in the inset. 
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FIG. 3. 1/Tr as a function of v. The lattice size L = 128 is used in the simulation. The dotted 
line with symbols correspond to the simulation results and the solid line shows the expected values 
from the definition (14). 
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FIG. 4. Tr as a function of L 2 . The dotted lines with symbols show the simulation results 
and the solid lines show the expected values from the definition (14). The four points in each case 
correspond to cases of L = 64, 128, 256, and 512. The simulation step in this figure is At = 0.2. 
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FIG. 5. Two simulation results for 1/Tr as a function of v, where the expected value of 1/Tr 
is also shown to guide the eyes. The lattice size and time steps used in the simulations are shown 
in the inset. 
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FIG. 6. Simulation and expected results of the velocity profiles, where L = 256, r = 2 {y = 2.5), 
t — ton = 600. The solid line shows the expected values from Eq. (13) with n = 1. The symbols 
show the simulation results. The type of the fluid and the value of t on are shown in the inset. 
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